Existence and nonexistence of an intrinsic tunneling time 
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Using a time operator, we define a tunneling time for a particle going through a barrier. This 
tunneling time is the average of the phase time introduced by other authors. In addition to the 
delay time caused by the resonances over the barrier, the present tunneling time is also affected by 
the branch point at the edge of the energy continuum. We find that when the particle energy is 
near the branch point, the tunneling time becomes strongly dependent on the width of the incoming 
wave packet, which implies that there is no intrinsic tunneling time. This effect is related to the 
quantum uncertainty in the particle's momentum. 
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I. INTRODUCTION 

The definition of tunneling time — the time it takes 
a particle to tunnel through a potential barrier — or 
even whether it can be defined or not, has been a much 
debated problem and is still a controversial one of fun- 
damental quantum mechanics E , H, [ p, |1 , d, H, 0, H, H, 
M M MM, M MMMm Mm- in this paper 
we address the question: Is there an intrinsic tunnel- 
ing time? We present a definition of tunneling time on 
the basis of a time operator canonically conjugate to the 
Hamiltonian [l^, [2l| . Our tunneling time consists of two 
contributions. The answer to the question is "Yes" for 
one contribution but "No" for the other contribution. 

The first contribution comes from the overlaps between 
the incoming wave and resonant states. It is basically a 
weighted sum over all resonance poles of the resonant 
lifetimes. We may say that the incoming wave splits into 
resonant channels of the tunneling barrier and spends 
the lifetime of each resonance before it tunnels out. For 
a particle represented by a spatially large wave packet, 
it is closely related to the phase time defined by Wigner, 
Smith, PoUak and Miller, and others 0, [H, [H, HI] . This 
is a dominant contribution to the tunneling time when 
the particle energy is near the resonance poles, It 
gives a tunneling time (as a function of the particle en- 
ergy) that is independent of the width of the incoming 
wave packet. In other words, it gives an intrinsic tunnel- 
ing time, which depends only on the resonance poles of 
the barrier. 

The second contribution appears when the particle en- 
ergy is near a branch point. In contrast to the first con- 
tribution, it is strongly dependent on the width of the 
incoming wave packet. This makes a universal definition 
of the tunneling time impossible near a branch point. 

In short, our main point is that when the energy of 
the incoming particle is near a resonance pole of the tun- 
neling barrier, an intrinsic tunneling time does exist, but 
when the energy is near the branch point there is no 
intrinsic tunneling time. This is understandable; while 



the resonance poles yield the Markovian dynamics (ex- 
ponential decay), the branch point yields non-Markovian 
dynamics (i.e., power-law decay) with no characteristic 
time or length scales, which cause deviations from ex- 
ponential decay for both long time scales ^24] and short 
time scales [l^HBl- In the following, we present general 
arguments to support our claim and present numerical 
results for a square-barrier model. 

II. TIME OPERATOR AND AGE 

Our argument starts with the time operator [l^, [2l| 



in imits with h = 1. Here H' is the part of the Hamil- 
tonian associated with a continuous spectrum, or the 
Hamiltonian excluding the bound states of the particle, 

-\E^)E^(E^\ (2) 

a -°° 

with denoting the dispersion relation of a mode a 
with wave number k of free propagation. In terms of the 
eigenstates of the Hamiltonian we have 




The time operator ([T|) satisfies the commutation rela- 
tions [t, H'] = i and [t, Hb] = 0, where Hi, is the part of H 
that includes the bound states. These commutation re- 
lations give e'^*te~'^* — t + t, so that the time evolution 
just adds time t to the time operator. This property al- 
lows us to interpret the time operator as giving the "age" 
of a state [13] . We define the average age of a normalized 
state \tp) as 

(i)^ = {i,m. (4) 
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Since the energy is bounded from below, the time oper- 
ator is not self-adjoint [H,!!!], i.e., {A\i\B) ^ {B\i\A)* . 
This means that the age of a given state at t = can 
be complex. However, we will consider the age difference 
between incoming and outgoing states of the particle. In 
momentum representation the states we consider will dif- 
fer only by a phase factor. As a result, the age difference 
will be real, despite the ages being complex. 

The difference in age between two states is given by 



(5) 



Keeping in mind an experimental scenario where the par- 
ticle has an average positive velocity (moving from left to 
right) and tunnels through a potential barrier, we define 
the initial state "01 as a state where the particle is known 
to be on the left side of the barrier, and the final state 
^2 as a state where the particle is known to be on the 
right side of the barrier, with both ^/^i and ^02 giving the 
same average velocity. Our postulate is that t^^,^-^ will 
then give an average of the time it takes the particle to 
move from the left of the barrier to the right. Note that 
■02 is not the time-evolved state ■01 (t) ■ If 02 were taken 
as 1^2 = i^iit), due to the relation e*^*fe~*^* = i + t, the 
age difference t^i(t),^i(o) would simply give t. 



III. GENERAL FORM OF THE AGE 
DIFFERENCE 

Let us consider a general one-dimensional system with 
a symmetric potential barrier centered at a: = 0. In posi- 
tion representation, the stationary eigenstates \Ek) (giv- 
ing the eigenvalue continuum) of the Hamiltonian have 
the form 

( r(A:)e^'=(^-°), X > a/2, 

{x\Ek) = < Bkix), -a/2 < a; < a/2, 

[ e**^^ -f i?(fc)e~*'=("=+"), X < -a/2, 

(6) 

where a is the width of the barrier, R is the reflection 
coefficient, T is the transmission coefficient, and Bk{x) 
is the wave function inside the barrier. 

We will use the symmetric {a = +) and anti-symmetric 
{a — — ) modes of the stationary states: 

^ ( ±e-''=^±i^±(fc)e*'=^, a:>a/2, 
{x\Ef) = - I Bk{x) ± Bk{-x), -a/2 < a; < a/2, 
^ y ^ikx _^ F±{k)e-'''^, X < -a/2. 

(7) 

The factor 1/2 accounts for the normalization of the 



states and the double-counting of positive and neg- 
ative k in Eq. The coefficient 

F± (fc) = {R{k) ± T{k))e-'*"' (8) 

is the scattering amplitude for the symmetric or antisym- 
metric outgoing waves, respectively. Since the incoming 
flux of e~*'^'l^l and the outgoing flux of e**^'^' should be 



equal for the stationary states, the scattering amplitude 
must have modulo 1, i.e., \F±{k)\'^ = 1 for real k. 

In position representation, we set the initial and final 

states as 



(x|^0i) = exp(iA:oa;)/v ^0 for — ig ^ fl/2 < a; < — a/2, 
(x|^02) = eyi^{ikQx) / ioi a/2 < x < Lq + a/2, (9) 

and for other x. These are truncated plane waves of 
width Lq. We have chosen truncated plane waves because 
we want to study the limiting situation when these states 
approach plane waves. This occurs when Lq ^ . 
Then both t/ji and ^02 approach plane waves with well- 
defined momentum /cq. We could use different functions 
such as Gaussians, but the functions above seem to be 
the simplest ones to consider. Hereafter we will refer to 
the truncated plane waves © as the wave packets. In 
momentum representation, we have 



(fci0i) - r(fc-fco), 

(fc|02) = e-*('=-'='')(^«+")/*(fc- fco), 
where f*{k — fcp) is the Fourier transform of {x\^Ji) 



(10) 



f*{k-ko) 



J{k-ko)a/2 



1 



pi(k-ka)Lo 



-i{k - ko) 



(11) 



In the end, we will take the limit — > cxd. This limit 
will turn out to be unique for some terms of the age dif- 
ference, but non-unique for other terms, because it will 
depend on other parameters such as the momentum of 
the incoming particle. We will judge the "intrinsicness" 
of the tunneling time by seeing whether it has a well- 
defined limit or not when the size of the wave packets 
tend to infinity, approaching plane waves with fixed mo- 
mentum. 

We will calculate the age difference ([5]) as i^2,'/'i = 
(V'2|t|V'2) - {ipi\i\ipi), where 



{ipj\i\iPj 



E 

Q = ± 



dk 

2^ 



. d 

''dEi 



imi^,)- (12) 



The age difference may be separated into the age differ- 
ence with no barrier, and a delay Ar due to the barrier. 



+ - -1- 



Ar. 



(13) 



The delay time At can be directly measurable, because 
it is just the difference between the average arrival time 
of the particle with the barrier present and the average 
arrival time with no barrier present. 

After some algebra (see Appendix |^ , we obtain 

tZ.^,=iL, + a)v-\ (14) 

where is the average inverse group velocity dk/dEk, 

dk 



f dk 

v{ko)-'^ / ^\f{k-ko)\' 



dEk' 



(15) 
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and Lq + a is the distance from the middle region of the 
initial state ipi to that of the final state ip2- Some more 
algebra (see Appendix|A| allows us to split the delay time 
into two parts: 



with 



At = Ata + Atb (16) 
,2 1 dF^W 



i,„w._!|:£|i/(,-,„,i^_ 



(fc) dEk 



(17) 



ArB(fco) = ^/^§(/(-fc + fco 



fi-k-ko)—fi-k + ko)) J^F^ik), 



(18) 



Q = ± 



which are real, because |F±(fc)p = 1, F±{k) = F^{-k), 
and /(k) = f*{—K) for real k and k. (In the calcula- 
tion we neglected terms of order 1/Lq). We will show 
in the next section that Ata is related to the motion of 
the particle inside the barrier, i.e., the tunneling process. 
On the other hand, Atb is related to the motion of the 
particle outside the barrier. We can see this because it 
only involves the reflection coefficient R, which is propor- 
tional to J2a=± -^a{k) in Eq. ^8]) : see also the discussion 
below Eq. §7^). 



IV. TUNNELING TIME 

In this section we will define the tunneling time. We 
will first relate Ata in Eq. (1171) to the phase time intro- 
duced by Wimer Q, Smith [2||, Pollak and Miller 
and others [1^| . This phase time is defined as 



Tph(fc)-5R(-*^lnr(fc) 



(19) 



for a particle with the wave number k. Writing T(k) — 
|T(fc)|e*''('=), we have Tph = 39 / dEk. To connect this 
to Ata, we write the amplitude as F±(k) = cxp(i0-|-(fc)), 
where 0±{k) is a real phase. Then 



Ata 



IE 



Q = ± 

Moreover, since 
1 



.(fc). (20) 



T = ^e'*^" (F+ - F_) = -e''"' (e*^+ - e'' 



(21) 
(22) 

we have e{k) = 7r/2 + {0+{k) + 6'_(fc) + 2fca)/2 and 

i ^ 6l„(fc) =e'(fc)-7r/2-fca, (23) 



tiinnel 



- average phase time 



age of \|/, 



age difference f<f^.if^ 
delay time Ax 



age of \|/2 



FIG. 1: Relation between the age difference, ttunnoi, ioutside, 
and the delay times Ata and Atb 



Therefore, from Eq. (I^Dl) we obtain 
Ata 



dk r) 



Using Eq. we finally obtain 
r°° dk 

^TA= ^\f{k-ko)\\i,{k)-av~\ (25) 

The fist term on the right-hand side represents a weighted 
average of the phase time, where the weight is the square 
modulus of the Fourier component of the incident wave 
packet. The second term av~^ is the time it would take 
the particle to move through a distance a (equal to the 
barrier width) by free propagation. The delay caused by 
the barrier is the difference between the tunneling time 
and the free-propagation time through the barrier region: 
Ata = ^tunnel — av^^. Hence we identify the timneling 
time with the average phase time as 



f'tuniicl 



av 



Ata 



see Fig. [T] 



/■°° dh 

/ f^l/(fc-fco)lV 
J-oo 27r 



EVALUATING THE INTEGRALS 



(26) 



In this section we will evaluate the integrals involved 
in the age difference, the delay time, and the tunneling 
time. We first note that in the momentum representation 
the time operator is given by 



t = i 



d 



d 



dEk k 



(27) 



which diverges at fc = 0. However, as we discuss now, this 
divergence is suppressed by the term {ipj \E^) in Eq. (fT2|) . 
which vanishes at fc = 0. In the k representation the 
eigenstates ([7]) are given by 



((Vi|fc)+F±(fc)(Vi|-fc)), 



{^2\Ei) = ±-{m-k)+F±{k){i,2\k)). (28) 
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When fc ^ we have F±{k) ^ — 1, because when k = 
the particle cannot tunnel through the barrier as long 
as the barrier has a positive height and positive width. 
The particle is perfectly reflected. Hence i?(0) = — 1 and 
r(0) = 0, which gives F±{0) = -1. (In Appendix [B] this 
is shown explicitly for a square-barrier potential). 

Thanks to this behavior of the scattering amplitude, 



we see that (V'jl-^'fc 



when k 



0, for J = 1,2. 



Moreover, the derivative d{E^\il}j) /dk is regular at fc = 
0. Thus the vanishing {ipj\E^) cancels the 1/k divergence 
at fc = coming from Eq. ([?7|) . This means that the point 
fc = is not a true singularity in the age difference. 

Since the integrand in the age difference is regular at 
fc = 0, we can replace 1/fc by its principal part without 
changing the integration. This means that we can replace 



1 
fc 



1 



1 



k + ie k — it 



(29) 



with e > real (infinitesimal). 

To evaluate the integrals in Eqs. (fT5|) and p6)) . we will 
also add an infinitesimal in the denominator of |/(fc — 
fco)P as (fc — fco)~^ {k — ko + ie)^^, which does not 
change the result because |/(fc— fco)P is regular at fc = fco. 
Similarly, we will add infinitesimals to the denominators 
of /(±fc ± fco) in the integral of Atb in Eq. P^ . 

Let us consider first Eq. ^T5\\ . or the average inverse 
velocity 



zvr 



dk 



dk 

2^ 



1 



l/(fc-fco)Py 

dfc If - e'(k-ka)Lo\ 



(30) 



m 



27r 



(fc - fco)2 



Expanding the absolute value squared and using Eq. (j29p . 
we obtain 



1 

m 



dk [1 



oj(fc-feo)Lol 



[1 



-i(fe-feo)Lol 



27r 



{k — kfj — iey 



(31) 



To evaluate the integral we will close the contour with an 
infinite semicircle either on the upper or the lower half fc- 
plane, depending on whether the integrand vanishes on 
the upper or lower infinite semicircle, respectively (see 
Fig. [2|). For the first term in brackets in the numerator, 
we will close the contour on the upper half plane, and 
for the second term in brackets, we will close the contour 
in the lower half plane. Evaluating the residues at the 
double pole fc = fco + ie and the poles fc = ±ie, we obtain 



m 
fco 
m 
fco fco 



1 



iknLo 



2Lo kl 
m sin(fcoio) 



2Ln 



fc2 



koLo 



(32) 





> 


Im(fe) 

. +ie K^i£ 


Re(jt) 




> 


; -IE 





FIG. 2: Integration contours and poles of the age difference. 



Next we evaluate the integral 
dk 



2tt 



\fik-ko)\'Tp^ik) 



(33) 



in Eq. ([26|1 . We follow the same procedure as for the 
calculation of v^^ outlined above. Taking the double 
pole at fc = fc'o + ie coming from |/(fc — fco)P and the 
poles at fc = ±ie coming from the derivative with respect 
to the energy in 



Tph 



d9{k) 



m 
2" 



1 



1 



dEk 2 \k + ie ' fc 
we obtain (with 9'{k) = de{k)/dk) 
dk. 



de{k) 
dk 



27r 



|/(fc-fco)|Vph(fc) 



(34) 



(35) 



= F^'(fco) 

fco 



m sin(fcoio) 



(fco) 



fco fco^o 
sin(fcoLo 



9'(0) + O(l/Lo) 



fc^Lo 



[fcTph(fc)]fe=o + O(l/Lo), 



where we neglected the residues at the poles of the 
phase time, which are also poles of the scattering ampli- 
tudes. As discussed in Appendix [Bl as long as the width 
and height of the barrier are non-zero, these poles give 
O(l/Lo) contributions, which we neglect for large Lq. 
(Note that the term involving sin(fcoLo) is non-negligible 
when fco ^ 1/^0-) Inserting Eq. (|35p into Eq. (|25|) we 
obtain. 



sin(fcoio) r, XI _i 
7-ph(fco) -aT~ ifc^ph(fc)Jfe=o - ■ 



k^Lo 



(36) 



For Atb we follow a similar procedure (see Appendix 
[C|). The result is 



fco' 



sin(fcoLo/2) 
fcoio/2 



(37) 
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The term mLQ/ko gives the time it takes the particle to 
travel the distance Lq /2 on each side of the barrier with 
a speed m/fco- 

The age difference is 



t 



1p2,1pl 



AO) 

V'2,V'l 



Ata + Atb 



Tph(fco) - ^^"fa^"'^"'' [fcT-ph(fc)]fe=o 



'-Lo 



This can be written as 



k^Lo 
sin(fco£o/2) 
fcoio/2 



^tunnel ~l~ ^outside; 



(38) 



(39) 



where 



itunnci = av ^ + Ata 

sin(A;oio) 



Tph(fco) 



and 



^'Outside 



Lqv-^ + Atb 



ko 



Lo 



sin(fcoLo/2) 
koLo/2 



(41) 



is the time that the particle spends outside the barrier. 



VI. BRANCH-POINT CONTRIBUTION 

In the previous section we evaluated the integrals in- 
volved in the age difference, by taking residues at the 
poles, including the poles k = ±je (with e 0). These 
poles are in fact associated with a branch point of the en- 
ergy. Indeed, the energy of the particle is Ek = fc^/(2m) 
outside the potential barrier. For the dispersion Ek oc fc^, 
the complex energy plane has two Riemann sheets Ek = 
k"^ /2m with Imk > and with Imk < 0. We have the 
branch cut on Ek > and the branch point at Ek = 0, 
or at fc = 0. For this reason we will call the residues at 
k = ±ie the branch-point contribution. 

The terms due to the branch point are the terms con- 
taining the sine function in Eqs. ((551) . 601), and (|1T|) . 
These terms vanish when kgLo 3> 1 but are non- 
negligible when koLo ^ 1. Since we are considering large 
wave packets, this means that the branch-point effect ap- 
pears when the momentum of the particle is close to zero 
(ko ^ 1/Lo). 

As /cq ^ and Lq oo the branch-point terms ap- 
proach either zero or infinity depending on the limiting 
order. In the limit limfc(,_,o limiQ^oo the branch-point 
terms vanish, but in the limit liuiLg^ac linifco_>o they di- 
verge. This extreme dependence on Lq means that the 
tunneling time has no characteristic scale near the branch 
point. It depends on the size of the incoming wave packet 



rather than any intrinsic time scale associated with the 
barrier. Because of this the tunneling time is not in- 
trinsic around the branch point. This gives a negative 
answer to the question of the existence of an intrinsic tun- 
neling time. As mentioned in Introduction, the branch- 
point contribution is associated with the non-Markovian 
dynamics (i.e., power-law decay) with no characteristic 
time or length scales 24, 25, 26]. 

We may understand the physical origin of the branch- 
point contribution as follows. The particle states that we 
are considering are, in position representation, truncated 
plane waves of large width Lq. In contrast, in momen- 
tum representation, these wave packets are high, narrow 
peaks centered at fco with a width or order 1/Lq. This 
width expresses the uncertainty in the momentum of the 
particle, which occur because the initial and final states 
are not plane waves but truncated plane waves. When 
ko approaches zero with ko < 1/ioj eventually be- 
comes smaller than the uncertainty range. As a result, it 
becomes increasingly likely for the particle to have neg- 
ative momenta. Negative momenta have the following 
effect on the age difference between two states: if two 
states ipi and '02 have negative momenta only, and "02 is 
located to the right of V"!, then -02 is actually "younger" 
than ipi. Therefore, negative momenta give negative con- 
tributions to the age difference between 02 and ipi. On 
the contrary, if the momenta are positive, then the age 
difference is positive. As a result, as fco ~* 0, the average 
age difference, including both negative and positive age 
differences, tends to zero. This can be verified by taking 
this limit in Eq. (|38l) . 

In short, the branch-point terms in the age difference, 
Eq. ([38]) . express a reverse flow of the particle due to 
momentum fluctuations rooted in the uncertainty princi- 
ple. The term involving sin(/co^o) represents the reverse 
flow through the potential barrier, while the term involv- 
ing sin^(fcoLo/2) represents the reverse flow outside the 
barrier. 

The branch-point terms give negative contributions to 
the age difference, and hence they decrease the average 
time it takes the particle to move from its initial state 0i 
to its flnal state ■02- In a sense, the branch-point causes 
the particle to speed up when it has a very small classical 
velocity. 



VII. RESONANCE CONTRIBUTION 

In this section we will isolate the contribution to the 
tunneling time coming from the resonances of the scat- 
tering amplitudes over the barrier. We will show that, in 
contrast to the branch-point contribution, this contribu- 
tion is an intrinsic function of Eq, being independent of 
the size of the wave packet. Moreover, we will show that 
the resonance contribution can be written as a weighted 
sum of lifetimes associated with each resonance of the 
scattering amplitudes. 
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The resonances appear in Ar^i; 

dk 



-y 



Taking the pole at k = ko, we obtain the residue 

d 



2 dE,„ 



In^^a(fco). 



(42) 



(43) 



We write the scattering amphtudes as functions of the 
energy Eq — /cg/(2TO) in the form 



Eo - {E 



En - E^ 



-G±{Eo), 



(44) 



where each F±{Eo) has two branches, since Eq — 
(±fco)^ / (2m). The product includes all the resonance 
poles E^ of the scattering amplitude. The complex con- 
jugate resonance (E^)* must appear in the numerator 
because F± has modulo 1 . The function G± is a remain- 
der factor with modulo 1 as well. 

Using Eq. HH) we obtain from Eq. US]) 



+ i 



■E 

Q = ± 

d 



E 



Eo - {Ef)* Eo - Ef 
\nG^{Eo) . 



(45) 



By writing the real and imaginary parts of the poles as 



(46) 



we obtain 



{Eo - Ef^y + (r^/2)2 



2EEt^ 

a=± 



(r,"/2)- 



(47) 



The first term in the right-hand side is the resonance- 
pole contribution to the tunneling time. This is the 
weighted average of the resonance lifetime with 
the Lorentzian weights between and 1. The factor 2 
in front of the summations takes account of the fact that 
the particle comes into and goes out of the barrier to 
tunnel, while the resonance lifetime is only the time it 
takes the particle to go out of the scattering potential. 



VIII. RESULTS FOR A SQUARE BARRIER 

In this section we will discuss the behavior of the tun- 
neling time and age difference as functions of the momen- 
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FIG. 3: The tunneling time ftunnci vs. the momentum of the 
particle fco, for a = 15, m = 1, 2mV = 1, Lo = 150 (solid), 
and Lq = 300 (dotted). 



tum of the particle fcg for £^ square barrier. We have 
The phase time is given [l^] by 



TO 1 ;nSinh(2Ka) + 2aKk'^(K'^ — fc^) , , 
Tnh(fc) = 7 o ^ (49) 

where l"^ = k"^ + k'^ = 2mV . Using Eqs. (gO]) and ((49)) 
we obtained the numerical plot of the tunneling time 
in Fig. [3] for two different widths of the incoming wave 
packet {Lo — 150 and Lo ~ 300), which are large com- 
pared to the width of the barrier (a — 15). As can be 
seen, the tunneling time has peaks near the resonance 
poles of the scattering amplitudes F± (around fco = 1.1). 
According to our theoretical analysis, in this region the 
tunneling time should be independent of Lq for large Lqj 
which is confirmed by the numerical plot. Thus itunnci is 
an intrinsic function of fco in the resonance region. 

On the other hand, near the branch point, i.e., near 
fco = 0, the tunneling time depends on the size Lo of 
the incoming wave packet. This means that there is no 
unique universal tunneling time. 

Figure [3] compares the age difference t^2,^i with the 
barrier present and the age difference t^^^ without the 
barrier. Note that around fco = 0.5 (inset), t^^^^-^ < 



t 



(0) 



. The particle is in a sense accelerated by the bar- 
rier. This is the Hartman effect [H, [13, [U, [13, [H . Close 

to fco = 0, i^s.j/'i is ^1^0 smaller than t^^l^^- This ef- 
fect is not related to the Hartman effect. It is due to 
the branch-point effect, that is, due to the reverse flow 
of the particle caused by momentum fluctuations, as we 
discussed at the end of Section IVII When fco is close to 
zero this reverse flow is enhanced by the reflection due to 
the barrier. This, on average, makes the particle arrive 
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FIG. 4: Mam plot: The age difFerence t^^.i'i (solid), and 
the age difference with no barrier i^^^ (dashed) vs. ko, for 
I/O = 150, a — 15, m = 1, and 2mV — 1. Inset: A zoom-in 
view around ko — 0.5. 



sooner when the barrier is present than when there is no 
barrier. 



Again, it will be important to detect the deviations in 
a tunneling experiment. 

Olkhovsky and Recami [2^ have argued that the do- 
main of the time operator should be restricted to func- 
tions of the energy E that vanish at i? = (the branch 
point). In this way the time operator becomes Hermi- 
tian and has real expectation values. Moreover, even 
if the domain includes functions that do not vanish at 
E — 0, they have proposed to use a bilinear time oper- 
ator, which, again gives real expectation values. In our 
approach, however, this issue is not very relevant because 
the age difFerence that we obtained is real to begin with. 
Using the bilinear time operator proposed by Olkhovsky 
and Recami [2§| instead of the operator in Eq. ([T]) will 
give the same age difference ([5]). 

As mentioned above, our tunneling time is the average 
phase time. The phase time has found some experimental 
support as a good measure of the tunneling time [ssj . 
(although it is not universally agreed that the phase time 
is the correct tunneling time). New experiments (and 
possibly more theoretical work) could confirm or negate 
our prediction. 
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IX. CONCLUDING REMARKS 

We have proposed a definition of tunneling time ob- 
tained from the change of expectation value of the time 
operator. We considered spatially large incoming and 
outgoing wave packets (truncated plane waves with av- 
erage momentum ko). Our tunneling time is the average 
phase time, averaged over the momemtum distribution of 
the incident particle. It reduces to the phase time Tph(fco) 
when the energy of the tunneling particle is far from the 
branch point of the energy continuum. 

However, near the branch point we obtain a devia- 
tion from the phase time Tph(fco); the deviation depends 
on the size of the incoming and outgoing wave packets. 
This deviation gives a non-intrinsic character to the tun- 
neling time. It may be interesting to see if a tunneling 
experiment of slow particles (with ko « l/^o) shows de- 
pendence of the tunneling time on the size of the wave 
packet as predicted here. 

Our calculations have centered on a symmetrical bar- 
rier, and symmetrical initial and final states. One could 
consider asymmetric configurations. Another possible 
extension of our work is to consider the case where the 
wave packets have positive-momentum components only. 
In this case we expect that there will still be branch-point 
effects, but they will take a form different from the one 
discussed in this paper. 

The definition of tunneling time we presented here is 
by no means the only possible definition. There are many 
other definitions, and one might wonder if our main pre- 
diction, the appearance of deviations from the phase time 
for slow particles, is not an artifact of our definition. 
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APPENDIX A: CALCULATION OF THE AGE 
DIFFERENCE 

Here we will show the main steps involved in the calcu- 
lation of the age difFerence between the final and initial 
states -02 and tpi. Defining 



(A^ivi) = r(fc-fco), 

(fclVa) = g*ik-ko)r{k-ko), 



where 



r{k-ko) = ^e 

V I^Q 



i{k-ko)a/2 



\ — gi(fe-feo)io 



and 



we have 



-i{k - ko) 



(Al) 
(A2) 



(A3) 



(A4) 



{^i\E^) = i,{f{k-ko)+F±{k)f{-k-ko)) 
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{^2\E^) = ±-{f{~k-ko)g{-k-kn) 
+ F±(k)f{k-k^)g(k-k^)) 

and (see Eq. (fT2l) ) 



(V-ilt^i 



Q = ± 



X [f{k-k^)—r{k-ko) 



(A5) 



(A6) 



as well as 



+F^{k)f{-k - k,)—F:{k)r{-k - fco) 

+F^{k)f{-k-ko)^J*{k~ko) 



+ f{k~ko)—F:{k)n-k-ko) 



i sr-^ f dk 



Q = ± 



X ( i^a(fc)/(fc - ko)g{k - k^)—F*^{k)g*{k - A:o)r (fc - fco) 



(A7) 



+ /(-fc - ko)g{^k - fco) /*(-fc - fco).g*(-fc - fco) 

oEk 

+ F^{k)f{k - fco)g(fc - ko)^r{-k - ko)g*{-k - fco) 
+ f{-k - ko)g{-k - ko)^Fl{k)r{k - ko)g*{k - fco) 



We first consider the first and second terms in the 
right-hand side of Eqs. (|A6p and (|A7p . Using the rela- 
tions |F±(fc)p = 1 and |5(fc)P ~ 1 for real fc, we find that 
they give the following contribution to the age difference: 



(A) 



■E 

a=± 



dk 

2^ 



(A8) 



d 



X [\f{k~kQ)\'g{k^kQ)—g*{k~ko) 



d 



+|/(-fc - fco)rg(-fc - fco)TT7^g*(-fc - fco) 



dEu 



For the terms involving — fc we change the variable of 
integration from fc to — fc. Noting that F^{~k) = F±{k) 
we obtain 



(A) 



dk rl 

^\fik-k,)\'g{k-ko)z^^g*{k-ko) 



Q — ± ^ 



c.c. 



(A9) 



where c.c. denotes complex conjugate. The first term 
in the right-hand side of Eq. (|A9[) does not involve the 



potential barrier. It therefore gives the age difference 
with no barrier, 



t 



(0) 



dk 3 

= I ^\fik - fco)Pff(fc - ko)^g^/{k - fco) 



(Lo + a) 
(Lo + a)v~^, 



^\nk-ko)\'^ 



where 



(AlO) 



(All) 



is the average inverse group velocity. On the other hand, 
the second term in the r.h.s. of Eq. (jA9[) may be written 
as 

a— ± 

because F± has modulo 1 and we have F|. = l/^±- 
When there is no potential barrier, we have F±(k) = ±1. 
Hence At/i vanishes when there is no barrier; it is a cor- 
rection to the age difference due to the barrier. 

For the third and fourth terms in Eqs. (|A6p and (IA7[) . 
we use the relation 



/(«)ff(K) = /*(«) 



(A13) 
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The contribution to the age difference from the third and 
fourth terms in Eqs. (|A6p and (|A7p is then 



(B) ^ 1 

Q = ± ■ 



d 



dk 



(A14) 



X (^F„(fc)r(fc-fco)^/(-fc-fco) 

+ ri-k-ko)^j:{k)f{k-ko) 
-F„(fc)/(-fc-fco)^r(fc-fco) 

-/(fc-M^^^o(fc)r(-fc-fco)) • (A15) 

Changing fc — > — fc in the second and fourth terms in the 
above expression and using F^{k) = F±{—k) as well as 
f*{k — fco) = f{—k + fcg), we obtain 



(B) • i f dk 

Q = ± 



X (^F„(fc)/(-fc + fco)^/(-fc-A:o) 

-F„(fc)/(-fc-fco)^/(-fc + fco)) . (A16) 

When there is no potential barrier, t\j,^ vanishes be- 
cause F±{k) — ±1 with no barrier. Therefore, ^^''^^ is 
another correction to the age difference coming from the 
barrier and we write it as Atb- The total age difference 
is 



(A17) 



APPENDIX B: SCATTERING AMPLITUDE FOR 
A SQUARE BARRIER 

In this appendix we will write down the scattering am- 
plitudes for a square barrier potential, and we will con- 
sider its limits when the momentum of the particle goes 
to zero and the width or height of the barrier go to zero, 
with the aim of establishing a range of validity of the age 
difference calculated in the text. 

The scattering amplitudes are obtained from their def- 
initions F±{k) — exp{—ika){R{k) ± T{k)) and the well- 
known expressions for the reflection coefhcient R{k) and 
the transmission coefficient T{k). The result is 



F±{k) - 



(1 ± e-'"')fc - ^(l T 6'"°)^^ 
(1 ± e"«")fc + i{l =F e-'=°)K' 



(Bl) 



where V is the height of the barrier, a is the width of the 
barrier, and k = \/2mV — k"^. 

We will consider first the limit a 0. It turns out that 
the limits fc — s- and a ^ of the scattering amplitude 
are not interchangeable. Indeed, when fc — > we have 



na = a^2mV . Thus, for any V > Q and a > 0, we have 
limfe^o F±{k) = — 1. However, when a ^ with fc > 
we have F±{k) ±1. We will discuss next how this 
affects the age difference. 

The integrands that appear in the age difference are 
peaked around fc = fco; therefore in the following discus- 
sion we will consider fc ^ fco . We will focus on the branch- 
point effect, which appears when fc ~ fco ^ I/Lq ^ 0. 
Moreover, we will consider the expression for the ampli- 
tudes when a is small, so that ko <C 1, or 



Then we have 



F+ik) 



F-{k) 



k + i{mV ~k'^/2)a 
k - i{mV - k^ /2)a 



k-2W2mV -k^/a 
k + 2W2mV -k'^/a 



(B2) 



(B3) 



(B4) 



These expressions show that for F_ the limits fc 
and a are interchangeable. For either limit we have 
—1. However, for F^ the limits are not inter- 
changeable, as mentioned earlier. When the limits are 
taken so that amV ^ fc ^ ^/Lq, then we have 



lim lim F+(k) 

a^Ofc^O 



-1, 



but when amV <ti k l/^O) then we have 
lim lim F+(k) = 1. 

fc^Oa->0 



(B5) 



(B6) 



This gives a discontinuity between the age difference with 
the barrier present (Eq. p8|) ) and the age difference with 
no barrier (Eq. (|14p) when we take the limit a — > in the 
former. This discontinuity appears because when we de- 
rived Eq. ([38]) we neglected the pole contributions coming 
from the scattering amplitudes, arguing that they gave 
0{1/Lq) corrections. This was fine as long as the width 
a of the barrier was finite. However, when a ^ 0, the 
poles of the scattering amplitudes give terms compara- 
ble to the term coming from the poles at k — ztie (the 
branch-point contributions) . 

Specifically, one can see that when a the scatter- 
ing amplitude F^{k) in Eq. (jB3p has a pole at fc = imVa. 
When a ^ this is essentially a pole at the branch 
point, giving a non-negligible residue. If we include this 
residue, the discontinuity mentioned above is removed. 
This brings the question of how large a has to be so that 
Eq. dSSl) is vahd. 

When a ^ 0, the residue of the pole at fc = imVa 
involves the term 



exp(ifcLo) ~ exTp{^mVaLo) 



(B7) 



coming from the incoming or outgoing wave-functions 
(see Eq. This term vanishes if 



a > l/{mVLo). 



(B8) 
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Therefore, Eq. §8l is valid only if Eq. (|B8j) holds. This 
condition is consistent with the condition we mentioned 
earlier above Eq. (jB5|l . 

Similar arguments apply when we take the limit V 
instead of a ^ 0. Equation p8)) is valid if Eq. (|B8p holds, 
i.e., if V ^ l/(maLo). 



cles. The integrals are then reduced to summations over 
the residues of the poles inside the contour. Since 



fi-k±ko 



1 



-i{k =F A:o) 



(C2) 



APPENDIX C: EXPLICIT FORM OF Ars 

In this appendix we will evaluate the term Ars in 
Eq. Writing the derivative d/dEk in terms of fc, 

we have (see Eq. (jA16l) ) 



--is/ 

a=± •' 



dk m 

2TT~k 

X (F„(fc)/(-fc + fco)^/(-fc-A:o) 

-F„(A;)/(-fc-fco)^/(-fc + fco)) . (CI) 

To evaluate the integral we will close the integration con- 
tour using either the upper or the lower infinite semicir- 



we will close the contour in the upper infinite semicir- 
cle. The functions /(— fc± fco) and their derivatives have 
no poles at A: = ±fco- The scattering amplitude F±{k) 
may have poles in the upper half-plane. However, any 
residues of these poles are, except for phase factors, in- 
dependent of Lq. Therefore, due to the I/^/Lq factor 
in Eq. (jC2p . the poles of the scattering amplitude give 
O{l/Lo) contributions, which we neglect. Finally, the 
1/k factor gives a branch-point contribution. Similarly 
to Eq. (j29p we interpret this factor as a principal part. 
Equation (|C1|) takes the explicit form 



Atb ^ Siko) - S{-ko), 



(C3) 



where 



Q = ± 



2 \k — ie k + ie 



or 



X e 



a = ± 
i{k+ko)a/2 



J(k-ka)a/2 



I _ gi(k-ka)Lo 
-i{k - fco) 

1 



d_ 
dk 



^i{k+ka)a/2 



^ _ ^i{k+kn)Ln 



2 \k ~ ie k + ie 



^i(k — ka)a/2 



i{k + fco) 

\ — gi(fe-feo)io 



_ ^i{k+ka)Lo 

—i{k + fco) 



1 



-iLo 



-i{k - fco) 



Jik+ko)Lo 



Taking the residue at fc = ie, we obtain 

1 - cos(fcoio) sin(fcoLo) 



2 k + koj -i{k + ko) 
I 



or, using cos(a;) = 1 — 2sin^(x/2) as well as Eq. 
arrive at 



klL. 



0-^0 



fc2 
ftp 



(C6) 



When fc = 0, the transmission coefficient is T = 
and the reflection coefhcient is R = —1. Therefore, the 
amplitudes F±(fc) = {R{k) ± T{k))e''"' become F±(0) = 
— 1 at fc = 0. Hence we have 



Atb = tJ^Q 
fco 



sin(fcoLo/2) 
fcoLo/2 



(C4) 



(C5) 



we 



Lov-\ (C8) 



Atb — —m 



1 ~ cos(fcoio) sin(fcoio) 



fc^Lo 



fc2 



, (C7) 
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